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^ ■ Abstract: The 2D quasi-geostrophic (QG) equation is a two dimensional 

■ model of the 3D incompressible Euler equations. When dissipation is in- 
S . eluded in the model then solutions always exist if the dissipation's wave 
>- ! number dependence is super-linear. Below this critical power the dissipation 

I appears to be insufficient. For instance, it is not known if the critical dis- 

■ sipative QG equation has global smooth solutions for arbitrary large initial 
■ — ' data. In this paper we prove existence and uniqueness of global classical 

solutions of the critical dissipative QG equation for initial data that have 
small L°° norm. The importance of an L°° smallness condition is due to the 
fact that L°° is a conserved norm for the non-dissipative QG equation and 
is non-increasing on all solutions of the dissipative QG., irrespective of size. 
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1 Introduction 

Do singularities develop in finite time in smooth solutions of unforced, incom- 
pressible 3D fluid equations? This challenging question remains yet unan- 
swered. Lower dimensional model equations have been proposed and studied 
([05 0)111 5 Hi) attempt to develop mathematical insight in this prob- 

lem. The 2D quasi-geostrophic (QG) equation is one of these models. The 
dissipative QG equation is 

39 

— + M- Ve + /€(-A)"e = 0, (1.1) 
ot 

where a G [0, 1], k > is the dissipative coefficient, and the 2D velocity field 
u = {ui,U2) is determined from by a stream function ip via the auxiliary 
relations 



Equation (LI) with a = 1/2 is the critical dissipative QG. Criticality means 
that the dissipation balances nonlinearity when one takes into account the 
conservation laws. 

In addition to its intrinsic mathematical interest, the equation ( |1 . 1| ) is 
relevant in the context of general quasi-geostrophic models of atmospheric 



and ocean fluid flow [10 



We are concerned here with global existence results for solutions of the 
initial-value problem (IVP) for equation ( p..lD wherein 

9{x,0) = 9oix) (1.3) 

is specified. We consider periodic boundary conditions with period box Q = 
[0,27r]^. Without loss of generality we may restrict the discussion to 9 that 
obey for all time (27r)~^ 9dx = 0. 

The issue of global existence for equation ( |1.1D is non-trivial (P],[P^]). 
If no smallness condition is imposed on the initial data, then the issue of 
global existence for arbitrary data is open. In this paper we show that if the 
L°°-norm of the initial data is small then ( |1 . 1| ) possesses a global solution in 
the critical case a = 1/2. The QG equations (dissipative or not) have global 
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weak solutions for arbitrary initial data (JTT|, see Appendix B). The L°° 
norm condition is significant in view of the fact that the QG equations have 
a maximum principle see Appendix A) that ensures that for all time 

and initial data the L°° norm is non-increasing in time. If the initial data 
is smooth enough then the solution of the critical dissipative QG is unique, 
smooth and decays in time. These results and their proofs are presented 
in section ^. When a > 1/2, the smallness assumption on the data is not 
needed for global existence of smooth solutions for equation (pTTI). This 
can be proved using the same ideas as for the critical case. The theory of 
global existence and regularity in the sub-critical {a > 1/2) case is thus in a 



satisfactory state. More details for the sub-critical case can be found in |14 



We establish now some of the notation. The Fourier transform of / is / 

[Zti)^ Jq 

A is used to denote the operator (— A)^, defined at the Fourier level by 

Af{k) = \k\f{k). 
The relation in (|1.2|) can be identified as 

u = {-d,,A-'e, d^,A-'e) = (-i?2^, RiO) or u{j) = I {j)^e{j), (1.4) 

where i = J = ^ for j G \ {0}, (ji, ^2)^ = (-J2, ji), and Ri and R2 

are Riesz transforms ([|12|)- The spaces are the familiar Sobolev spaces 
of functions having s derivatives in L^. 



2 Global existence 

The initial value problem for the critical dissipative QG equation is 

' Ot + u-ve + nAe = 0, {x,t) eQ X [0,00), 

< u = {ui,U2) = {-R2O, RiO), (x,t) G X [0,00), (2.1) 
^ e{x,0) = Ooix), X eQ, 
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where k > is a constant. 

In this section we assume that the norm of the initial data is smalL 
We estabhsh that the I VP (|2.1|) has a global bounded solution in H^. If the 
initial data is smoother (i?^) then the solution's norm in is non-increasing 
in time. The solution becomes real analytic at positive time and and decays 
exponentially. For the sub-critical case, no smallness assumption is necessary 
and these results hold for arbitrary data. 

We start with an apriori estimate. 

Theorem 2.1 There exists a constant Coo such that for any Oq G fl 
with 

\\0o\\l^ < Coo/t, (2.2) 
the classical solution 9 of the IVP (\2.1\) satisfies 



m-,t)\\H2 < \Mh2 (2.3) 

for all t > 0. 

Proof. Multiplying the first equation in ( |2 . 1| ) by and integrating by 
parts, we obtain 



1 d 
2'dt 



J \Ae\Ux + K J \{-A)-^efdx = - J A^e{u-ve)dx. 

Further integration by parts gives 

J A^e {u ■ ve)dx = 2 jvu- {yiye))Ae + j{Au- ve)Ae 



By Holder's inequality. 



A^9 {u ■ Vd)dx 



< C 



\Vu\\L4Ad\\i, + \\Au\\L4Vd\\L4Ad\\Li 

The Riesz transforms are bounded in W spaces, so 

\\Au\\l3 < C||A^||i3, llVnlliS < C\\Ve\\L3. 
The Gagliardo-Nirenberg inequalities 

llv^lUs < c\\e\\i^\\{-A)U\\l,, \\Ae\\L. < c||^^||i^||(-A)t^|||„ 
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follow from classical ones (0) by complex interpolation (see also[Q). Using 
them we obtain 



< c\\d\\Loo\\{-A)^^e 



|2 

Il2- 



J {-Af9{u ■ V9)dx 
Collecting the above estimates, we have 

J \A0\Ux + K J \{-A)-ie\^dx <Coo\\0\\L^\\{-A)-*e\\l2. 



1 d 

2dt 



It was proved in ([0) (see the Appendix A) that 9 satisfies the maximum 
principle 

||^(-,i)||L- < ll^olU- foralH>0. 

Taking Cqo = (C'oo)"^ the bound ( |2.3| ) then follows from the smallness condi- 
tion (|2.2| ). This completes the proof of the theorem. 



Theorem 2.2 There exists a constant c^o (the same as in Theorem 
above) so that for any 9q G with ||^o||l°° ^ Coo/t the IVP l \2.1\ ) has a 
unique global solution 9 satisfying 

\\0{;t)y2 < \\9o\\h^ 

for any t > 0. 

Proof Let 9 be the unique local solution on [0, Tq] with Tq depending on 
1 1 6*0 1 1 only (standard techniques can be applied to show that a unique local 
solution exists and depends continuously on initial data in H^). By Theorem 
pTTI , 9 satisfies 

\\0{;t)\\H^ < \\9o\\h^ 

for any t G [0,To]. Therefore the local solution can be extended uniquely to 
[0, 2Tq] and the global solution is obtained by repeating this procedure. 

Theorem 2.3 Assume that the initial data 9q G satsifies the bound 
||^o|U°° < CoqK. Then the solution of the IVP ( \2.1\ ) decays exponentially 

II 6' II ^2 < exp (—ct) 11^0 11^2 

for all t > 0. Here c = 2{k, — cJ^^||^^o||l°°); one? Coo is the same as in the 
previous theorems. 
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The proof is a trivial consequence of the obvious Poincare inequahty 

/ |(-A)t^|2rfx> / |(-A)^|2dx. 
Jn Jn 

Theorem 2.4 Assume that the initial data 60 G satsifies the bound 
||^o||l°° < CoqH- Then there exists a time to > such that the solution of 
the I VP ( \2.1\ ) is real analytic for t > to. More precisely, there exists an ex- 
tension of 9, 0{z,t) that is an analytic function in the time- expanding strip 
c C2, 

T,t = {z = X -\- iy; x E Q, \y\ < -^{t — to)} 



and obeys the inequality 



\0{z,t)\<'^ 



uniformly for 2; G . 

Proof. Let us consider to to be the first time when 

Y{t)= E mt)\ 

jez2\{o} 

becomes smaller than k/A: 

Y{to) < k/4. 

In view of the preceding theorem and the elementary inequality 

Y{t)<c\\e{t)\\H2 

the existence of to is guaranteed. Consider the function 

{t-to)^\j\\ 



j6Z2\{0} 
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and the function 

{t-to)K\j\' 



4t)= E \j\mt)\exp^ 

J6Z2\{0} 
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Formally 

| + .W(f-yW)<o 
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holds. This imphes that the function y{t) is non- increasing for t > to. The 
more rigorous proof requires one to take only a finite sum, and introduce an 
artificial power to avoid differentiating the modulus at zero: 

ynAt)= E l^(j,t)r+^exp 

ieZ2\{0},|j|<n 

This is now a differentiable function in time. One differentiates, and obtains 

d K 

-gyn,e{t) + -Zn,e{t) < ^,y{t)z{t) 

where 

Zn,e{t)= E |j||^(j,t)|l+^exp 

ieZ2\{0},|i|<n 

and 7e = (1 + e)||6'o||^2, and thus limj^oTe = 1- One integrates from t = to 
to t and passes to the limit n —>■ oo, e ^ 0. One obtains then 

y{t) + 77 / z{s)ds < 7 + / y{s)z{s)ds 

Then, because y{to) < f it follows that the set {t > to\ y{s) < I, Ws,to < 
s <t} equals [to, cc). The completely rigorous proof requires a regularization 
of the equation so that z{t) is guaranteed to be finite. After t = to one may 
use Galerkin approximations for this purpose. The conclusion is that, for 
t >to one has 

iez\{o} { ^ ) ^ 

The analytic extension is 

e{z,t)= J2 e''-nj,t) 

jez2\{o} 

and the uniform convergence and bound in follow. This completes the 
proof. The idea of using time dependent exponential weights was introduced 
in (0). 

Initial data in are sufficient for a global existence results to hold. 



{t-to)K\j\ 



{t-to)K\j\ 
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Theorem 2.5 There exists a constant d^Q such that for any 9q G and 

||^o||l°° < doot^ there exists a weak solution of the QG equation satisfying 

\\e{-Mm < WOoWm 

for any t > 0. 

Proof It was proved in ([|lT|) that the IVP (|2l] ) with 6q G has a global 
weak solution 0{-,t) G L"^ (see Appendix B). The weak solutions are 
constructed using a Galerkin approximation. The weak solutions can be 
constructed by solving approximate equations 

dtO + U5-Ve + kM = 0, (2.4) 

where 

Us = ks * u = ks * (-R"*"^) 
and ks the periodic Poisson kernel in 2D given at the Fourier level by 

ksiO = e-'l«l 

^ G Z^. The approximations have gobal smooth solutions for positive time, 
uniform bounds in L°°{dt; L'^{dx)) and converge weakly to solutions of the QG 
equation. In addition, and in contrast with Galerkin approximations, these 
approximations have monotonic non-increasing norms of u. Multiplying 
the first equation in (|2.4| ) by and integrating by parts, we obtain 

1 d 
2~dt 



J \ve\'^dx + kJ \Ale\'^dx < J \{ve) ■ vus ■ {ve)\dx. 



By a similar argument as in the proof of Theorem |2.1| , the term on the right 
hand side can be bounded as follows. 

J \{v9)-vus-{V9)\dx < ciiv^iiis ||v«5||l3 < c\\ve\\ls < cii^iUooIiaI^ii^^. 

In addition to the boundedness of Riesz transforms in spaces and an appro- 
priate Gagliardo-Nirenberg inequality, we use here the fact that convolution 
with the Poisson kernel does not increase norms. Therefore, 



ld_ 
2di 



J \ve\'^dx + K J \A^e\'^dx < D^\\e\\Loo\\Ak\\l2. 



8 



Using c/oo = (-Doo) ^ and the maximum principle we deduce from ||6'o||loo < 
(ioo and the inequahty above that 

\\9{-,t)\\Hi <\\0o\\h^ and J \A'k\'^{x, T)dx dr < oo (2.5) 

for any t > 0. The solutions obtained thus are relatively strong and the 
equation holds in time integral form in L^. 

For the sub-critical case a > | the maximum principle allows one to get 
a sublinear bound of the nonlinearity in terms of the dissipation and, conse- 
quently, the global existence result holds without any smallness assumption 
on ll^lli^oo. 

Theorem 2.6 Let a > and 6q G . Then there exists a unique global 
solution 6 solving the IVP ( \2.1\ ). The solution is real analytic for positive 
time and decays exponentially to zero. 



3 Appendix A 



In UTTI it was shown that a solution to for | < a < 1 satisfies the 



following maximum principle 

||^(-,^)IUp < II^oIIlp forl<p<oo foralH>0. 

Below we give a short description of the proof in the case a = |. 
Define 9 = kg * 9 where ks{^) = e"'''^'. The Poisson kernel kg is positive 
and has integral equal to one. 9 satisfies 

as 



and 

ds 

Integrating with respect to s we obtain 

rs2 



+ pj \9\P-^9A9dx = 0. 
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Using the properties of the Poisson Kernel kg 

ll^si * 0\\Ip - \\ks^-s^ * {ks, * 6)\\Ip > 0. 

therefore 

J \e\P-^eAedx = /im,^o+ / \0\^'^OAedx > 0. 
Hence, returning to the evolution equation for 9 



d\\9 



IP 

I LP 



dt 

4 Appendix B 



-P J M^'^dAddx < 0. 



It was proved in ([0) that the QG equations, dissipative or not, have global 
weak solutions in L^. We present here a brief description of the main reason 
for this fact. We consider QG as an infinite system of ODEs: 



^ e{i, t) + t) = bi{e, 6) (4.6) 



dt 
with 

He, 6)= ^{j^-k)eij,t)eik,t) (4.7) 

j+k=i u\ 

Note that because of the perpendicularity 

{j^-k) = {j^-l)=-{k^-l) 
and because of symmetry considerations in the sum we can write 

bi{e,e)= y: iikOij,t)mt) (4.8) 

j+k=l 



with 



\k 






\m 





Note the inequality 

< 



2(max|j|; \k\) 
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Consider now the weak norm 

||^|U= sup |^(i)| 
jez2\{o} 

The nonhnearity can be written as 

iGZ2\{0} 

The main observation is that the nonhnearity has a weak continuity property. 
Let 9 1 and 62 be in L^. There exists a constant C so that 

\\A-\B{9,,9,)-B{92,92))\U< 

C\\9, - ^slU (1 + log (1 + 11^1 - ^2||;')) (||^i||l2 + IMl^) 

This is the key ingredient in the existence of weak solutions. The method of 
proof employs a Galerkin approximation. The weak continuity guarantees the 
fact that the weak limit of the approximations solves the weak formulation 
of the equation. 
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